Adiabatic pumping in a double-dot Cooper-pair beam splitter 
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We study adiabatic pumping through a double quantum dot coupled to normal and supercon- 
ducting leads. For this purpose a perturbation expansion in the tunnel coupling between the dots 
and the normal leads is performed and processes underlying the pumping current are discussed. 
Features of crossed Andreev reflection are investigated in the average pumped charge and related to 
local Andreev reflection in a single quantum dot. In order to distinguish Cooper-pair splitting from 
quasi-particle pumping, we compare the properties of Cooper-pair pumping with single-electron 
pumping in a system with only normal leads. The dependence on the average dot level and the cou- 
pling asymmetry turn out to be the main distinguishing features. This is contrasted with the linear 
conductance for which it is more difficult to distinguish single-particle from Cooper-pair transport. 

PACS numbers: 74.45.-|-c,73.23.Hk,72.10.Bg 
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Charge transport through interfaces between super- 
conductors and normal conductors takes place by dif- 
ferent processes. Above the energy gap of the super- 
conductor's density of states, mainly single electrons are 
transferred, while subgap transport is sustained by An- 
dreev reflection (AR)ii In an AR process, an electron 
in the normal conductor that impinges on the interface 
is retroreflected as a hole while a Cooper-pair is trans- 
ferred into the superconductor. In junctions with more 
than one normal conductor also crossed Andreev reflec- 
tion (CAR) may occur, that is, the two electrons forming 
the Cooper-pair stem from different normal conductors 
(or tunnel into different normal contacts in the opposite 
transport process). This nonlocal transport mechanism 
has been extensively studied both theoreticallji^ and ex- 
perimentallyi^ 

In recent years advancements in nanofabrication have 
made it possible to contact quantum dots (QDs) with su- 
perconducting lcads4r— Such QD-superconductor devices 
are of great relevance, because they enable the investi- 
gation of the interplay between superconducting corre- 
lations and Coulomb repulsion in nonequilibrium situa- 
tions. Andreev reflection^"— as well as crossed Andreev 
reflectioniir— through quantum dots have been the fo- 
cus of many theoretical works. Recently, CAR through 
QDs has also been observed in experiment 1^1^ The setup 
consists of a superconducting lead tunnel coupled to two 
parallel quantum dots realized in an InAs semiconducting 
nanowiro^ and a carbon nanotube^ respectively. Each of 
the two quantum dots is additionally coupled to separate 
normal reservoirs (see Fig. [1]). As a result the depen- 
dence of the current in one arm of the beam splitter on 
the parameters of the other arm indicates the occurrence 
of CAR. 

In the examples mentioned above, a bias voltage is ap- 




FIG. 1: (Color online) NDSDN setup: two quantum dots 
coupled to the same superconductor and each dot coupled to 
a normal conductor. 



plied to generate dc transport. In this paper we consider 
a different transport mechanism: adiabatic pumping. The 
principle of pumping is to transport electrons in the ab- 
sence of a bias voltage by varying certain system param- 
eters periodically in time. Pumping is therefore a mech- 
anism converting an ac into a dc signal, which has been 
experimentally realized in different systems j^^ In the adi- 
abatic regime the pumping period is large compared to 
other characteristic time scales of the system. It was 
showni^ii^ that adiabatic pumping reveals features which 
are not visible in stationary transport. Here, the main 
motivation of our work is to use adiabatic pumping in or- 
der to investigate features of CAR. We therefore consider 
a system consisting of two quantum dots coupled both to 
normal and superconducting leads as shown in Fig. [T] In 
the experiments performed so far, the CAR and AR sig- 
nals coexist. This happens even though strong Coulomb 
interaction within each dot tends to suppress AR, there- 
fore enhancing the visibility of CAR. Adiabatic pumping 
requires two out-of-phase time-dependent parameters in 



order to obtain a finite dc current. Choosing gates ap- 
plied to the two dots, belonging to two different arms of 
the beam splitter, as pumping parameters, only transport 
mechanisms relying on nonlocal correlations between the 
two arms contribute to the pumped charge. Therefore, 
such a type of pumping cycle has the advantage with re- 
spect to biased transport that it singles out CAR, while 
local effects do not yield any finite dc current. 
Theoretically, the dynamical scattering approach pro- 
vides a general framework for pumping as long as the 
Coulomb interaction is weaki^i In noninteracting sys- 
tems, the influence of the superconduc ting proximity 
effect on pumping was studied in Refs. [l^ However, 
Coulomb interaction cannot be neglected in the setup 
considered here. In recent years much effort has been 
put on the treatment of pumping through strongly inter- 
acting systems such as quantum dotsji^ While pumping 
through a single quantum dot with a superconducting 
lead was studied in the limit of zero temperature and 
infinitely strong Coulomb interaction;^ in this paper, 
we are interested in systems in which Coulomb interac- 
tion can be arbitrary and where coupling to the leads 
is weak. To this purpose we use an adiabatic extension 
of a generalized master equation approach] ^^'^^i^^'^^ In 
the stationary limit the generalized master equation ap- 
proach^^ has been applied to hybrid quantum dot sys- 
tems before^iiiii^ 

The motivation of this paper is to identify and un- 
derstand CAR in adiabatic pumping. To this purpose 
we investigate two quantum dots, with infinite intra- 
dot Coulomb repulsion, tunnel coupled to the same su- 
perconductor and each of them to a normal conductor 
(NDSDN) (see Fig. [1]). Pumping is realized by apply- 
ing time-dependent potentials, namely one to each of the 
quantum dots, via gates with a phase-difference in the 
driving. This gives us the possibility to identify unique 
features of crossed Andreev reflection in adiabatic pump- 
ing which rely on the nonlocality of the effect and can - 
as we show by a comparison with a setup with the super- 
conductor replaced by a normal lead (NDNDN) - not be 
reproduced by other parasitic nonlocal effects mediated 
by quasiparticles. 

However, the complexity of this setup makes it difficult 
to obtain compact analytic formulae. Therefore, we addi- 
tionally consider a quantum dot with Zeeman-split levels, 
tunnel coupled to a ferromagnetic and a superconducting 
lead (FDS). In hybrid systems containing ferromagnets, 
superconductors, and quantum dots the influence of spin 
asymmetry on Andreev reflections has been investigated 
beforei^^— 1^ In the present work our motivation of con- 
sidering the FDS setup is to get a better understanding 
of the transport processes in the NDSDN system because 
we can relate the CAR in the NDSDN setup to AR in the 
FDS setup. The Zeeman splitting and the polarization 
in the FDS setup corresponds to a difference of the two 
dot levels and an asymmetry of the coupling to the two 
normal conductors of the NDSDN system, respectively. 
From a theoretical point of view the main difference be- 



tween the two setups is the existence of triplet states 
in the NDSDN system. Experimentally, although hy- 
brid systems containing ferromagnets and superconduc- 
tors are realizable^^ the time dependence of the trans- 
port channels through the dot are easier to control in the 
NDSDN setup. 

The paper is structured as follows. In Sec.|ll]we present 
the models of the considered setups. The technique used 
to compute the pumping current into the superconduc- 
tor is described in Sec. IIIII The results divide in three 
different parts. In Sees . II V !^ and HV B I the results for lo- 
cal Andreev reflections and crossed Andreev reflections, 
respectively, are given. The features of CAR and single- 
particle transport are compared in Sec. I IV CI Finally, 
conclusions are drawn in Sec. IVl 



II. MODEL 

The systems we consider are generally described by a 
Hamiltonian for a hybrid system composed by multiple 
quantum dots tunnel coupled to both normal and super- 
conducting leads. Each individual dot, labeled by the in- 
dex r, is described by the Anderson-impurity model with 
an onsite interaction C/intra and the level energy Era-- The 
interaction between electrons in different dots is charac- 
terized by the inter-dot repulsion C/intcr- The quantum 
dots are described by the Hamiltonian 
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where n^a = dl^dra is the number operator for electrons 
in the dot r with spin a and dra (dla-) is the correspond- 
ing annihilation (creation) operator. Here we explicitly 
introduce the time-dependence of the dot levels, which 
is used to realize the pumping cycles. The leads are de- 
scribed by the Hamiltonian 
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where the different reservoirs are identified by the label 
ij. The operator Crfka {cl^ka) annihilates (creates) an elec- 
tron with momentum k and spin a in lead 77. The second 
term is only present for the superconducting leads and it 
is simply the attractive potential of the mean-field BCS 
Hamiltonian. Without loss of generality the pair poten- 
tial A can be chosen to be real, because wc consider only 
one superconductor. Finally, the dots are coupled to the 
different leads by means of the tunneling Hamiltonian 



-'^tiinn 



'tunn — /_^ '-W^r^ku 
rjrka 



^>ir'^nkrr"'rcr + h.i 



(3) 



Both the tunnel matrix elements and the density of states 
of the leads p^ are chosen to be energy independent 



in the window relevant for transport. Tunnel- coupling 
strengths are then defined as T^^r,a = 2TT\t^^r\'^pr,,a- No- 
tice that no inter-dot tunneling is included in the model. 
Finally, the total Hamiltonian for this type of hybrid sys- 
tem can be written as H = H^ot + -fftunn + J2-n-^v We 
set in the following h = I. 



A. Double-dot device 

The main focus of this paper is on the parallel double- 
dot device shown in Fig. [I] that is ideal for studying 
Cooper-pair splitting. It is composed of two quantum 
dots which are tunnel-coupled to different normal con- 
ductors but the same superconducting lead. We will re- 
fer to it as to the NDSDN system, where N indicates a 
normal lead, S a superconducting lead and D a quantum 
dot. The Hamiltonian of the NDSDN system is obtained 
from the general Hamiltonian of the previous subsection 
by having r G {L,R}, rj £ Nl,Nii,S and Tn^ = ^Nl.l, 
r^Vn = ^Nr.r, ^s,r as spin-indcpcndcnt tunnel-coupling 
strengths. With this we define F^r = FjVj^ -|-F7v„. For the 
double-dot system we assume the dots' levels to be spin 
degenerate, that is e^^ = £^4, = £r, the Coulomb repul- 
sion within one dot to be infinite C/mtra — ^ 00, and a finite 
inter-dot interaction [/inter = U. The limit C/intra — > 00 
excludes the possibility of double occupation of the same 
dot and, therefore, only CAR and no local AR appears. 
As independent pumping parameters we choose the two 
spin-degenerate dot levels, {eL,£_R}, which can be var- 
ied by means of gate voltages. This system will be con- 
trasted to the system with the lead S in its normal state, 
which is referred to as NDNDN and in which we take 
T]eNL,NR,Ne. 



B. Single-dot device 



neling can be obtained by integrating out the supercon- 
ducting degrees of freedomi^J^i^^ Here we will discuss 
the resulting effective Hamiltonian only for the NDSDN 
system. The one for the FDS system is completely analo- 
gous. The effective Hamiltonian in the limit f/intra ~^ 00 
readsi^ 



HeB = 2_/ ^rnr<y + C^ 2_/ '"■LariRa' 
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with Ts = V^'siXsR. being the effec- 

tive coupling. The eigenstates are |x) G 

{|+),h),|a,0),|0,a),|T_i),|To),|ri)} , where 

|o-, 0) (|0, cr)) corresponds to the left (right) dot be- 
ing singly occupied with spin a and the right (left) 
dot being empty. The triplet states are jT-i) = |i,i), 
l^o) = (|;,t) + |t,;))/V2 and |Ti) = |t,t>. The 
tunnel-coupling to the superconductor leads to eigen- 
states that are coherent superpositions of the state 
with both dots empty |0, 0) and the singlet state 

\s) = i\u)-n,i))/V2, 
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where S = Sl + sr. + U is the detuning between the 
empty state and the singlet and 2eA = yj ^"^ + 2F| is the 
energy splitting between the |-|-) and |— ) states. The 
corresponding eigenenergics arc E± = | ± £^, ii'(cr,o) = 
£l, -E'(o,cr) = £fl, and Et_^ = Etq — Et^ = £l + £_r + C/. 
In the FDS setup the eigenenergics and eigenstates are 
the same except that L and R are replaced by f and J,, 
respectively, the triplet state does not exist, the singlet 
state 15") is replaced by a double occupation \d) ~ did] |0) 
of the dot, and 2eA 



V^ 



r|. 



In order to identify the processes relevant for pumping, 
we consider a single-level quantum dot tunnel coupled to 
a ferromagnct and a superconductor (FDS), which having 
a smaller Hilbert space allows for a simpler analysis. The 
Hamiltonian of the single-dot system is obtained from the 
general Hamiltonian considering only one dot (we consis- 
tently drop the index r) and two leads: rj & F,S. The fer- 
romagnct is described by the Stoner model which induces 
F-|- 7^ F|. The nonvanishing tunnel-coupling strengths 
are: Tp and Ts- The pumping cycle in this case is re- 
alized by varying independently the two spin-split levels 
£-|-,£4.. This can be done by means of a time-dependent 
gate voltage and magnetic field. 



C. Large- A limit 

In the A — > 00 limit quasi-particle transport in the 
superconducting lead is suppressed and an effective de- 
scription of the dot that takes into account Andreev tun- 



III. METHOD 

A. Generalized master equation 

The system, which is described by the Hamiltonian 
given in the previous section, can be subdivided into two 
different subsystems, the (proximized) quantum dots and 
the normal conducting leads. Since we are not interested 
in the dynamics of the leads' degrees of freedom, we can 
trace them out. This leads to an effective description of 
the quantum dots in terms of the reduced density ma- 
trix prcd- The elements of this reduced density matrix 
are denoted by p^i ^ (xilProd IX2), where xi and X2 
are states of the dots. The diagonal elements p-^ 
give the probability to find the dots in state x- We in 
troduce the vector n — (p^^, ...,p^^^, ...,p^'., ...)""", with 
i ^ j, where the first m components are the diagonal ele- 
ments of the reduced density matrix of an m-dimensional 
Hilbert space followed by the off-diagonal elements. The 
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dynamics of the reduced density matrix is governed by a 
generalized master equation (in matrix notation) 



dt 



by means of Fermi's Golden Rule. Solely for the ones con- 
necting off-diagonal elements this is not sufficient and one 
has to resort to a diagrammatic method which has been 
developed in Refs. [2j. In general, ofFdiagonal elements of 
7r(i) = -iE(t)7r(i) + j dt'W{t,t')TT{t') . (6) ^j^g reduced density matrix, p^,, enter Eqs. and ©. 

However, we assume weak coupling to the normal con- 
ductors Fjv ^ ksT, ea, where for the FDS as well as the 
NDSDN setup the offdiagonal elements of the reduced 
density matrix are decoupled from the diagonal onesi^ 
As we are interested in the diagonal elements, needed for 
the computation of the current, we can therefore disre- 
gard the offdiagonal ones. Solely in the NDNDN setup 
the dynamics of the offdiagonal elements P(q'J) and pL'g( 



The matrix elements Wyy,,{t,t') of the kernel describe 
transitions from an initial state at time t' described by 
Py,, to a final state at time t described by p^,. In the 
systems which we consider, consisting of a single dot or 
two dots coupled in parallel, the matrix elements of E(t) 
are given by E^^^,„{t) = S^^^S^^^n, {E^{t)-E^,{t)). 

We study transport relying on the periodic variation 
of a set of pumping parameters {Xi(t)}. Assuming the 
parameter modulation to be slow, that is the pumping 
frequency fl to be small compared to all other energies 
of the system, we can perform an adiabatic expansion 
of Eq. dni) following the lines of Ref. [l5|. Within the 
adiabatic expansion with respect to a reference time t, 
the reduced density matrix is written as the sum of an 
instantaneous contribution and its adiabatic correction, 
7r(t) — > -Kf + TTj . The instantaneous contribution re- 
sults from freezing all parameters to their value at time 
t and yields a contribution in zeroth order in il/F^v, in- 
dicated by the superscript (i). The fact that the actual 
state of the system always slightly lags behind the pa- 
rameter modulation is captured in the adiabatic term of 
first order in i7/Fjv, indicated by (a). 
On top of the adiabatic expansion we perform a system- 
atic expansion in the weak tunnel-coupling strengths be- 
tween normal conductor and leads, F^r < ksT, of the ker- 
nel and the reduced density matrix, taking into account 
tunneling processes up to first order in F^r. Orders in the 
perturbation expansion in the tunneling coupling are de- 
noted by numbers in the superscript. The instantaneous 
contribution to the reduced density matrix is determined 
by 
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together with the normalization condition mr 
with n = (1, ...1,0, ..., 0), that is, the first m compo- 
nents of n are 1 and the other components are 0. The 
Laplace transform of the Kernel at zero frequency, with 
all parameters frozen to the time t is given by Wj = 
lim f (it'e~^(*~* ^Wl (i — t'), where here we consider 

only the first order in T^ (if not specified otherwise) . The 
adiabatic correction to the reduced density matrix turns 
out to have a contribution in minus first order in Tpf^ 
which due to the adiabaticity condition fl/T^ ^ 1 is not 
divergent and it is given by 



dt 
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is coupled with the dynamics of the occupation proba- 
bilities. In the NDNDN system, where also ofFdiagonal 
elements of the reduced density matrix contribute, we as- 
sume Ae = El — Efl ~ Fat and E and Wj ' have to be 
of the same order in the small parameter Fjv — Ae^^^ 
In a similar way, one can write rate equations for the 
expectation value of the current into lead ij. The instan- 
taneous contribution to the current is 



/«(t) = enW;''''V 
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(9) 



which we consider in first order in the tunnel cou- 
pling, only. The current rates W^' take into account 
the number of electrons transferred to lead rj^ From 
Eq. (O, we derive the conductance, which is given by 
G = (d/(*'i)/dy)|y=o, with V being the bias voltage. 
The instantaneous current vanishes exactly in the ab- 
sence of an applied bias. The adiabatic correction to the 
current is then the dominant one and it is given by 



/i;'°)(i) = enWj''(''^V 
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We are in the following interested in the charge trans- 
ferred into lead 77 per cycle of the parameter variation. 
This is found by integrating the current over one period 
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with nTr'"' ^' = 0. The rates W^ between diagonal 

elements of the reduced density matrix can be obtained Jn^'^' dt6Xi{t)-^6X2{t) 



In the following we consequently drop the index ?/ if 
the pumped charge corresponds to the superconductor, 
Qxi,X2 = Qxi x^- Two time-dependent parameters are 
necessary to create a nonvanishing pumped charge; we in- 
dicate the parameter choice in the subscript. The pump- 
ing parameters can be written as Xi (t) = Xi + SXi (t) , 
where Xi is the mean value and SXi{t) the oscillating 
component. We concentrate on the limit of weak pump- 
ing, that is, the oscillating component is small compared 
to the tunnel coupling SXi{t) ^ Fat. Therefore, we 
only account for terms up to bilinear order in SXi(t) 
and the pumped charge is proportional to Axi,X2 = 
2iT/n 



IV. RESULTS 



Local Andreev reflection 



Using the effective Hamiltonian and performing the 
perturbation expansion as presented in the previous sec- 
tion we calculate the pumped charge in lowest order in 
Tn or Tp, respectively. Close to the dot levels being at 
resonance, the lowest order processes are the dominant 
ones and cotunncling processes can safely be neglected. 
Before tackling the more complicated problem of CAR, 
we will first study the FDS system, in order to under- 
stand the features of local AR in adiabatic pumping and 
to identify the different transport processes occurring in 
this simple setup. For this setup, we also examine the in- 
fluence of cotunncling processes on the pumped charge 
far from resonance (Coulomb-blockade regime), which 
are important when the interaction U becomes much 
larger than the temperature. In Sec. IIVBI we discuss 
how adiabatic pumping provides the possibility to study 
CAR. To this end, we finally compare the NDSDN setup 
with the NDNDN setup. 



In this section we consider adiabatic pumping through 
the FDS setup. We choose the dot-level positions for 
electrons with different spins e^{t) and ei{t) to be the 
pumping parameters. Such a situation can be realized 
by a time-dependent gate voltage and a time-dependent 
magnetic field, the latter introducing a time-dependent 
Zeeman splitting. This choice of pumping parameters is 
convenient here as it allows for a direct comparison with 
a double dot in the absence of a magnetic field, in which 
gate voltages applied to the two dots are independently 
modulated. Pumping is possible whenever the polariza- 
tion of the leads or the average level splitting Ae = e^—ei 
are nonvanishing. To get a better understanding of the 
transport properties we first focus on two different limits: 
a vanishing polarization {p = 0) and a vanishing average 
level splitting {e^ = e_i). We start with the case of a 
vanishing polarization and a finite level splitting. 

For the pumped charge we find 
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with f'{x) — -^f{x) being the derivative of the Fermi 
function. We made use of the approximation f{E^ — 
£t) « f{E+ - ej « and /(e^ - E+) « /(e; - E+) « 1, 
which is justified for Fg > ksT. Equation ([T2|) shows 
that the pumped charge vanishes for an average Zeeman 
splitting equal to zero, that is e-j- = ej.. 

In Fig.[2l^a), we show the pumped charge Qei-e,: with- 
out the approximation on the Fermi functions used to 
write Eq. ((T2|) . as function of the average value of the 
mean dot level e = (e^ -I- £i)/2. The pumped charge 
exhibits a three-peak structure. The two external peaks 
are observed when the dot is in resonance with the nor- 
mal lead, that is when the addition energy for a single 
electron equals the chemical potential. This is realized 
for E- — Eo- = 0. Since we consider Zeeman splitting Ae 
being larger than ksT (with e^ being the level with the 
lower energy), only the resonance E- — ej^ = is accessi- 
ble due to Coulomb blockade. The other Andreev bound 
state, with energy E+ , is only accessible in the high-bias 
or high temperature regime. The two resonances associ- 
ated to the condition £'_ — £4^ = are at the two posi- 
tions, £^axa± « ^(-U±[{U + \Ae\f - F|]i/2), and are 
enhanced for an increased average Zeeman splitting. 

The central peak appears when the dot is in resonance 
with the superconductor, that is the average dot level is 
emax,2 ~ —U/2 which is realized for S < T^. In this case 
the dot undergoes fast oscillations between the empty 
and doubly-occupied state due to coherent Cooper-pair 
transfer. In particular these oscillations are much faster 



than tunneling events of single particles between the nor- 
mal conductor and the dot. However, transport requires 
exchange of charge both with the normal and the su- 
perconducting leads. Therefore, increasing the Coulomb 
repulsion U leads to an overall suppression of the pumped 
charge. The three peaks are not suppressed in the same 
manner. The side peaks are suppressed by the factor 

[1 + (C/ + e^ + £_^)Vr|]~', appearing in Eq. dH]). In- 
stead the central peak is suppressed by the combination 
of Fermi functions in Eq. p^ . 

We now focus on the limit of a vanishing average level 
splitting (A£ = 0) and a finite polarization. The pumped 
charge is then given by 



Qet,e,(A£ = 0) 



approximating the Fermi functions as done above. We 
find that the pumped charge is an odd function of S, 
therefore vanishing at the electron- hole symmetric point. 
The full result for the pumped charge at zero average 
detuning, A£ = 0, is shown in Fig. [H^b). As shown in 
Eq. ((T3| . the pumped charge vanishes at £ = ~U/2. 
However, we find a peak- trough structure, that is, the 
maximum contribution to the pumped charge appears in 
two peaks, close to £ w —U/2, with opposite sign. As 
argued above this relies on fast Cooper-pair oscillation. 
The amplitude of the pumped charge is much larger than 
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finite polarization, p ^ 0, we show it in Fig. ^c) as a 
function of tlie average dot level. The shape is a com- 
bination of the two structures shown in Figs. [2ja) and 
[2Ib). We find that the effect for the finite polarization 
dominates. Therefore, the peaks around £ sa —U/2 with 
opposite sign are the main feature to identify the proxi- 
mization of the dot. 

When lowering the temperature, the height of the 
peak-trough structure increases with inverse tempera- 
ture, that is, it becomes more and more pronounced. 
This result can, however, only be trusted as long as 
temperature is still large enough such that all charge 
states are thermally occupied. In the Coulomb-blockade 
regime, U » ksT and 6 = ef + ei + U ~ k^T, when the 
sequential tunneling rates to reach an empty or doubly- 
occupied dot are exponentially small, higher-order pro- 
cesses such as cotunneling need to be taken into account. 
To compare with the result presented in Eq. (|13p , we an- 
alyze the pumped charge in the cotunneling regime. For 
this, we first of all note that Eq. ([8]) looses its validity 
in the Coulomb-blockade regime, since the rates VF^*'^) 
get exponentially suppressed, while - in contrast to situ- 
ations where the magnetic field is constaniii^ - the time- 
derivative of the instantaneous occupation probabilities, 
^P(^ , of single occupation with spin a do not. The 
time-evolution of the probabilities of single occupation 
is then governed by spin-fiip processes in second order in 
the tunneling, W^.J, entering Eqs. ^ and ([TU]) together 
with adiabatic corrections to the probability in minus sec- 
ond order in F, pfjj . However, since U ^ UbT and 
5 PS fceT results in an exponential suppression of ^Pt±\ 

also the elements p^ j. are suppressed and will not enter 
the current in the Coulomb blockade regime. For the cal- 
culation of the cotunneling rates we follow the procedure 
introduced in Refs. [28| for metallic islands and applied 
for single-level quantum dots, for example, in Ref. |29| . In 
contrast to Eq. pop in the cotunneling regime the current 
is then 
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FIG. 2: (Color online) Pumped charge Q = Qe^.e, in units 



of Qo = 



eA^ 



(ksTy 



as a function of the average dot level e. The 



parameters in all figures are Fg — AksT and U — lOksT. 



which is nonvanishing due to p)°2 — 






= -Pil 



and 



Due to charge conservation, 

Qe^.si — '^Qe e • ^^^ charge pumped into the super- 
conductor is found as 



ief,e_i 



in Fig. Ufa) and strongly depends on the polarization of 
the leads: the stronger the polarization the larger the am- 
plitude. Furthermore, the pumped charge, in the vicinity 
of the electron-hole symmetric point is not suppressed 
by the strong Coulomb repulsion. We will address this, 
when discussing the cotunneling regime. 

Instead of giving the lengthy expression of the pumped 
charge for a finite average level splitting, Ae ^ 0, and a 



Zen^kBTA,^.,j:lp (1 - p2) 5 (r| + ^2 _ [;2) 



(1 +p2)7r2(fcBT)2F| + ^(1 - p2) (r| -f ^2 _ Jj2y 



(15) 



where we used Ae/C/ ^ 1. The qualitative behavior of 
the pumped charge in the cotunneling regime strongly 
differs from the sequential tunneling regime. For strong 
Coulomb interaction, in the cotunneling regime transport 



is suppressed with 1/U^ . To find a possible explanation 
for this suppression we focus on the transport processes 
during one pumping cycle. Consider the following pro- 
cess where a net transport is obtained in the cotunneling 
as well as in the sequential tunneling regimes: An elec- 
tron tunnels from the ferromagnet onto a singly occupied 
dot. The dot is then, for example, in state |— ). To obtain 
a net transport another electron has to tunnel from the 
ferromagnet onto the quantum dot bringing it back into 
single occupation which is possible due to Cooper-pair 
oscillations. A comparison of the system's time scales 
for the two regimes might shed light on the origin of 
the suppression of the pumped charge. In the sequential 
tunneling regime the time between two single-electron 
transport processes scales with l/F^v. In the cotunnel- 
ing regime the intermediate state can only be virtually 
occupied due to energy conservation and hence the time 
between two tunneling events scales with 1/U . In the 
considered limit of large U ^ fcsT and small Tn ^ fc^T, 
Cooper-pair oscillations are fast compared to the time be- 
tween two tunneling events in the sequential but slow in 
the cotunneling regime. This gives an interpretation of 
the suppression of the pumped charge in the cotunneling 
regime. 
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B. Crossed Andreev reflection 

We now consider a system made out of two quantum 
dots each coupled to one normal conducting lead. The 
two QDs are then coupled to each other via a common 
superconducting lead (see Fig. [T]). We take the pair 
potential in the superconducting lead to be the largest 
energy scale (A — > oo), such that single-particle trans- 
port between superconductor and QDs is suppressed. 
Furthermore, we take the intra-dot Coulomb repulsion 
(t^intra -^ oo) to be large excluding double occupation of 
each of the single dots, as discussed in Sec. Ill CI In this 
regime, only nonlocal effects enable transport between 
the superconductor and the dots, that is a Cooper-pair 
has to be split into two electrons occupying different dots 
or electrons from different dots enter the superconductor 
to form a Cooper-pair. 

We now calculate the charge pumped through the system 
due to the periodic modulation of the dot levels ei,{t) 
and s-[{{t), which can be achieved by two time-dependent 
gate voltages. We are interested in the charge, Q^^.en, 
pumped into the superconducting lead, which due to 
charge conservation and to the fact that only CAR is 
allowed is twice the charge pumped out of each normal 
lead. 

In Figs. [2Ja) and [3{b) we show Q^^^cr ^^ a- function 
of £ for different values of Ae and for different cou- 
pling asymmetries with the normal conducting leads, 
A = (Tnl ~ ^N!i)/^N, respectively. Features appear 
at the resonance condition with the normal and super- 
conducting leads, that are equivalent to the one in the 
FDS case with Zeeman splitting replaced by the differ- 



FIG. 3: (Color online) Pumped charge Q = Qe^.s^ in units 



of Qo = 



eA^ 



(kBTy 



as a function of the average dot level e. The 



parameters are Ts ~ SksT and U = lOfesT. 
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FIG. 4: (Color online) Linear conductance as a function of 
the average dot level e. for different coupling asymmetries A. 
The parameters are Ts = SksT, Tn = ksT, U = WksT, 
and Ae = ksT. 



ence of the energy levels of the left and right dots and the 
polarization p replaced by coupling asymmetry A. If the 
couplings to the normal leads are symmetric, A = 0, the 
charge as a function of the average mean dot level po- 
sition e, shows three peaks similarly to the FDS case. 
In this respect, CAR exhibits similar features to AR 



through the single dot. The main difference between the 
two is the asymmetry in the heights of the external peaks 
which can be attributed to the triplet blockade discussed 
in Ref. |l3J. Since the proximization by the superconduc- 
tor solely causes a coupling between the empty and the 
singlet state, Cooper-pair tunneling is blocked whenever 
the dot is in the triplet state. In the FDS setup the 
symmetry of the two external peaks can be related to 
particle-hole symmetry which is broken by this triplet 
blockade in the NDSDN structure. 

As in the FDS with finite polarization, also in 
the NDSDN the scenario changes completely in the 
asymmetric-coupling case (A 7^ 0). In this case the peak 
at e = — C//2 is replaced by a large peak-trough struc- 
ture. Interestingly, this feature dominates the external 
peaks which are barely visible in Fig.[3l^b). The position 
of the maxima and minima of this feature are exchanged 
when reversing the coupling asymmetry (A — ^ ~A). 

However, in the linear conductance, the coupling asym- 
metry does not introduce any new feature, as shown in 
Fig. 31 where for different coupling asymmetries only the 
weight of the three peaks is influenced and not their po- 
larity. Furthermore, the central peak is strongly sup- 
pressed. That means that the characteristic features of 
CAR in adiabatic pumping are not present in the linear 
conductance. As we will see in the next section these fea- 
tures are fundamental to distinguish single-particle trans- 
port from CAR. 



C. Single-particle transport 

A finite pumped charge can be obtained by varying in 
time the properties of the two spatially-separated dots 
exclusively by nonlocal correlations. CAR has such a 
nonlocal character. However, there may be other nonlo- 
cal effects that can produce a finite pumped charge and, 
thus, mask the signal from CAR. In order to distinguish 
CAR from other nonlocal transport processes, we investi- 
gate single-particle transport in a NDNDN setup, where 
the superconductor in the NDSDN setup is replaced by a 
normal conductor. While in the NDSDN setup the non- 
locality arises from CAR, in the NDNDN setup pumping 
is possible due to the formation of a coherent superposi- 
tion of states with one electron either in the left or the 
right dot. This superposition is generated by the tun- 
nel coupling to the common normal lead. In contrast to 
the NDSDN setup, the coherent superposition is strongly 
suppressed if the difference of the two dot levels is large 
compared to temperature (|Ae| 3> fcsT). 

Furthermore, in the NDSDN setup pumping cannot 
lead to an average charge transfer from the left into 
the right normal lead (and vice versa) because transport 
through the superconductor always involves CAR in the 
infinite- A limit. Instead, in the NDNDN setup, charge 
can also be transferred from the left lead Nl to the right 
lead Nr. Therefore an asymmetry of transport into lead 
Nl and into lead iVjj is one possible indication for single- 
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(Color online) Pumped charge Q = Q^,e^ in units 



eA, 



(ksT) 



^ as a function of the average dot level e 



FIG. 5: 

of Qo = 

for different coupling asymmetries A. The other parameters 
are Tiv.l = 0.4/csr, Fat^ j, = O-^ksT, Fn = O.lfcsT, U = 
WksT, and Ae = kBT. 



particle transport. 

The motivation of this work is the identification of 
CAR with respect to quasi-particle transport in form of 
an easily detectable signature in the pumped charge. We 
find this to be the peak-trough structure at e = — C//2 
that appears in the presence of a coupling asymmetry. In 
single-particle transport, modeled by the NDNDN setup, 
this feature is completely absent and only the peaks at 
the normal resonances appear (see Fig. [5]). These nor- 
mal resonances also have opposite signs, which cannot 
be reversed by changing the coupling asymmetry. Also 
in the very special situation of a symmetric coupling 
(A = 0) CAR can be distinguished from single-particle 
transport by the presence of the peak at e = —U/2. 
Therefore, an experimental study of the pumped charge 
in the double-dot system as a function of A as well as its 
behavior around e « —U/2 can clearly distinguish CAR 
from quasiparticle transport. 



V. CONCLUSIONS 

We have investigated adiabatic pumping through two 
quantum dots tunnel coupled to the same superconduc- 
tor and additionally coupled to different normal conduc- 
tors. For an infinite intra-dot Coulomb repulsion in this 
setup pumping relics on CAR. In order to understand 
the underlying transport processes we mapped the setup 
to the simpler setup of a quantum dot tunnel coupled 
to a ferromagnet and a superconductor where only AR 
appears. We found that most of the features of pump- 
ing including CAR are also present in pumping with lo- 
cal AR. The main difference are asymmetries due to the 
presence of the triplet state. To distinguish CAR from 
single-electron tunneling, which does not appear in our 



model but might be relevant in experiments, we com- 
pare transport through the double-dot setup containing 
a superconductor with a setup where the superconductor 
is replaced by a normal conductor. The dependence on 
the average dot-level position and the dependence on the 
coupling asymmetry A turn out to be the main distin- 



guishing features. 
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